function f = loss_LR(lambda_n)
%input: column vector lambda_(nth iter)
%output: calculated cost value
global Hpos Hneg;
f = sum(log(1+(exp(-(Hpos*lambda_n)))')) + sum(log(1+(exp((Hneg*lambda_n)))'));




